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Abstract 

In this short note we study the fourth order consequences of the stochastic sta- 



bility property for mean field spin glass models introduced in [AC]. We show 
that due to a remarkable cancellation mechanism it reduces to the well known 
second order version as predicted by Parisi's replica- equivalent ansatz. 
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1 Introduction 



In recent times a considerable amount of effort has been devoted to clarify the statistical 
mechanics structure of those models that describe disordered systems like the spin glass. 
A common feature of the considered class of models is the presence of two time scales, 
with different order of magnitude, governing the microscopic dynamics which reflects 
on a precise ensemble prescription in equilibrium statistical mechanics commonly called 
quenched ensemble. 

Among the different approaches to the problem the so called mean field considers 
systems in which every particle interacts with the remaining ones with the same random 
strength (see JSg, P?V| , [[TAIJ and @). 

In |[AC|| a stability property for the mean field models was investigated and its 
meaning in terms of continuity (in the temperature) of the thermodynamic functions 
was clarified. In particular some of the consequences of such a stability were proved to 
be captured by an infinite family of constraint for the overlap distributions. The result, 
in full agreement with the Parisi solution, provided an infinite family of overlap poly- 
nomials with zero average. Similar results were obtained in Q and later generalized 
GG). 



in 



The stability property, has been then investigated in [[FMPP1|| and [|FMPP2|| and 
cleverly used to determine a relation between the off-equilibrium dynamics and the 
static properties. Among other things such a relation allows an experimental (and 
numerical) measurement of the equilibrium order parameter for spin glass models and 
answers positively, at least in the mean field case, the important paramount question 
of whether the equilibrium properties of a glassy system have any relation at all with 
the off-equilibrium properties, the only ones physically accessible. 

It is then of crucial importance to have a clear understanding of the stability prop- 
erty. The study can be approached at least from two different points of view. One, in 
the framework of the GREM models, is the study of the random measures over fields 
of points on the line with the property of invariance under a reshuffling dynamics (see 



0, H] and 

The second, that we approach in this paper, deals with the characterization of the 
stability property in terms of the overlaps distribution. 

The paper is organized as follow: section 2 introduces the model and states the 
main result in proposition (2.3). Section 3 contains its proof. 
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2 Definitions and Results 



The SK spin glass model ||SK|| JMPVj has spin variables <7j = ±1, i = 1, . . . , N, inter- 
acting via the Hamiltonian 



H N (J,a) 



(2.1) 



l<i<j<N 



where the Jij are N ^ N ~^ independent normal Gaussian variables. 

Using the symbol Av to denote the Gaussian integration it is possible to prove [|GT 
that the quenched free energy density defined by 



-Pf N (P) = jjAvQ°gZ(p,J)) 



(2.2) 



converges (increasing monotonically) in the N — > oo limit. 

If we rescale an independent sample of the energy H by a factor \/N we define the 

K N (J',a) = -1 J l^3 , (2-3) 

l<i<j<N 

where J-j are i.i.d. unit Gaussian independent from J i; j of which turns out to be a 
Gaussian family indexed by a whose covariance is the positive definite matrix: 

-I 2 



2/ i\ 

q {a, a) 



1 \- , 



(2.4) 



It is extremely useful to introduce the quenched measure over an arbitrary number 
of copies of the system. For such a purpose we may consider r copies of the same 
system: each system is given the Hamiltonian where I ranges from 1 to r. To each 
copy we associate the Random-Gibbs-Bolztmann state given by the usual formula: 

j2 - e -P H{l) i°) 



< - > 



a) 



^ e -^H(i)(,r) 



(2.5) 



<< - >> 



< - > 



(0 



(2.6) 



The full measure E(-) is then defined as 



E(-) = Av{®\ =1 < - > w ) 
The matrix Q, whose elements are defined by: 

Qi,m = q 2 (a il \cr irn) ) 



(2.7) 



(2- 
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is then a symmetric random matrix with respect to the measure E(-) because its 
distribution turns out to be invariant under the action of the permutation group over 
the r copies. For instance for r = 3 



(2.9) 



1 


Ql,2 


Ql,3 




l 


Q2,3 


Qs,i 


Q3,2 


1 



and 



E(QwQw) = Av{[Z((3,E)]- 3 Yl gVV (2) )gVV (3) ; 

'}. (2.10) 



CT (l) i(7 (2) )(T (3) 



We define as in [|AC| the Random- Gibbs-Boltzmann deformed state as 

< - exp XK > 

< ~ >XK := ~Z \ 7? ~ 5 

< exp XK > 

<->ak := ®!<->$r, (2.11) 

and 

E XK (-) := Av{<£ - >ax) • (2.12) 

It may be shown that the Wick integration formula for Gaussian systems imply that 
all the thermal observables of the SK model (i.e the derivatives of the free energy 
with respect to (3) are indeed moments of the Random matrix Q with respect to its 
permutation invariant measure E. 

In [|AQ] it was shown the stochastic stability property: 



Proposition 2.1 If, over a certain temperature range, the states E^ N '^(—) are uni- 
formly continuous in beta, as N — > oo, then the limiting measures are stable 
under deformations in XK, in the sense: 

E(-) = Exk(-). (2.13) 



Always in [ ACfl it was derived as a consequence of the stochastic stability that an 



infinite family of overlap polynomial has zero mean. In fact it was proved that: 
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Proposition 2.2 At finite values of N and for every overlap monomial M 

' -E N (M) = N—E NiXK (M)\ x=0 = NE N (AM), (2.14) 



(3d(3 JVV ' d\ 2 ' 

where AM is a suitable overlap polynomial. In particular, if the quenched state is 
locally differentiable in f3 uniformly in N , the expectation values of quantities of the 
form AM vanish at the rate 0(1/N) or faster i.e. in the limit 

d 1 



d\ 2 



E XK (M)\ X=0 = E(AM) = 0. (2.15) 



Remark: the family of polynomials AM whose mean is zero plays in the mean field 
spin glass models a role which is reminiscent of the Hermite polynomials in classical 
Gaussian theories; in the same way the identities that can be found in \GQ\j allow a 
computation of the moments in terms of sum of products of lower order ones and are, 
in parallel, reminiscent of the Wick formula for the Gaussian case. 

In this short note we want to investigate if stochastic stability when considered 
for higher order derivatives in A lead to new results with respect to its second order 
consequences as one may a priori expect or if it is rather already captured by the second 
order as conjectured by Parisi in his replica equivalent ansatz (a generalization of the 
replica symmetry breaking ansatz |F|]). For this purposes we consider here the fourth 
order derivative with respect to A (the odd-order derivatives are zero for parity, the 
general even-order case will be considered in |pK||) and we show the following results: 



Proposition 2.3 For every overlap monomial M 

d 1 



<9A 4 



E XK (M)\ X=0 = 3E(A 2 M). (2.16) 



As a corollary of ( 3.25Q and (|2.16| ) and considering the fact that AM is a polynomial 



in the overlap algebra i.e. AM = c a ^iM a the vanishing of the average of AM for 
every M implies also the vanishing of the average of A 2 M: 

— E XK {M)\ X=0 = 3E{A 2 M) =3^c a , AM £(AM a ) = , (2.17) 



where the last equality holds due to the (|2.15|). In conclusion the consequences of 



stochastic stability at the fourth order are included in the consequences at the second 
order as predicted by the replica- equivalent Parisi ansatz (JPJ. An easy but important 
corollary of the previous proposition is the following 
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Proposition 2.4 In / |AC| / it was proved that 



(M en(m) = n ' 




<9A 4 



E N , XK (M)\ X=0 , 



(2.18) 



due to the ( 2.11 ) we have 



d 



) 



2 



(38(3 



E N (M) = 3N 2 E N (A 2 M) . 



(2.19) 



If, over a certain temperature range, the states En(-) are uniformly differentiable in 
beta (Ci), as N — > oo ; then over the same range (up to, at most, a countable number 
of points) they are uniformly (C2) and the expectation values of quantities of the form 
A 2 M vanish at the rate 0(1/N 2 ) or faster. 

3 Proof of proposition 2.3 

In the sequel we will have to deal with derivatives of functions involving Gaussian 
variables with covariance Q. The functions we will consider will be products of Gaus- 
sian and the Boltzmann weights. A typical quantity could be Qf 2 Q2,3- Computing 
the A-derivatives of E{Q\ 2 Q2,z)\ will amount to consider repeated truncated expec- 
tations. It may be useful then to introduce a graphical representation in which a 
monomial of the previous form is identified with a graph whose vertexes are the replica 
indices {1, 2, 3, ...r} and the edges correspond to the overlaps Qi tTn . Such a graph will 
be indicated by the symbol (1,2) 2 (2,3). The treatment of unpaired Gaussian vari- 
ables will be denoted by free legs and denoted by symbols like (1), (2) etc. It is 
important to stress that for each monomial in the overlap variables the invariance un- 
der permutation of (E, Q) will reflect on invariance under permutation labeling over 
graphs: (1,2)(2, 3) = (2, 1) (1, 3) = (2, 15) (15, 3). Moreover we will only consider the 
(I, I) = 1 case, which covers all the p spin models, and Derrida's REM and GREM 
structures. On a mathematical point of view it is then convenient the introduction 
of the notion of generalized graphs which is an object containing both edges and legs: 
(1,2) 2 (3), (1, 2)(3, 4) 2 (1). The permutation symmetry over indices will imply for in- 





(3.20) 



v€V(G) 



5 V G = 5^ + <fH , 



(3.21) 
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5^G = (v)G , (3.22) 

Ex: ^ +) (1,2) = (2)(1,2). 

S^G = -{v')G (3.23) 

where v' is any index not belonging to G. 

Ex: 6^ (1,3) = -(2)(1,3), ^"'(l, 3) = (4)(2)(1,3), 

We introduce then the Wick contraction map C which acts only on the legs part of 
the graph according to the Wick formula for Gaussian 

Let / be set of indices with even cardinality, then 

c n« = e um- (3.24) 

idl pairings v,v' £l 

Ex: C[(l)(2)(3)(4)] = 3(1,2) 

Proposition 3.5 For every generalized graph G 

d 4 



E Xh (G)\ x=0 = E(C5 4 G), (3.25) 



<9A 4 

where CS A G is a polynomial in the matrix entries. 

The proof of the above proposition is straightforward. The operation S is the graphical 
counterpart of the usual derivative with respect to the parameter A in the Boltzmann 
weight (where it appears in XK). Such a derivative produces a truncated correlation 
expressed in the rule ( |3.2i| ). The ( |3.2U| ) is nothing but the Leibniz rule for derivative 
of products. Each differentiation with respect to A produces a sum of monomials, each 
containing an added centered Gaussian variable (K) of zero mean. The contraction is 
nothing but the Wick rule for normalized Gaussian families. 

We will get the result ( |3.25D as a direct consequence of the following combinatorial 



theorem: 

Proposition 3.6 Defining A := CS 2 , 

CS 4 = 3A 2 . (3.26) 

Remark: it is interesting to note that the left hand side of the last equation contains a 
priori a number of addends which is 3-16 = 48 (the 3 coming from the Wick contraction 
of a fourth order monomial and the 16 coming from 2 4 terms of (5 + + 5~) A ). The 
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right hand side instead contains 4 • 4 = 16 terms from the square of A. Although 
the Wick contraction doesn't conserve the number of edges nor the number of vertexes 
(see formulas 3.44 to 3.47 in the next section) the presence of alternating signs in 
the definition of 5 together with the invariance of permutation of graphs labeling will 
produce a delicate cancellation mechanism responsible for the clean equality ( \3.6j ). 

Proof. Let show, for the sake of simplicity, the theorem for the moment E(Q 12 )- 
The generalization will be discussed at the end. 
As in [ AC | we have 

^A(l,2) = (1,2) 2 - 4(1,2)(2,3) + 3(1,2)(3,4) (3.27) 
Applying the A to the previous formula addends we obtain 

iA[(l,2) 2 ] = (1,2) 3 - 4(1,2) 2 (1,3) + 6(1,2)(3,4) , (3.28) 

^A[(l,2)(2,3)] = + 2(1,2) 2 (2,3) - 9(1, 2)(2, 3)(3, 4) + (1, 2)(2, 3)(3, 1) 

- 3(1,2)(1,3)(1,4) + 6(1,2)(1,3)(4,5) (3.29) 



^A[(l,2)(3,4)] = + 2(1,2) 2 (3,4) + 4(1, 2)(2, 3)(3, 4) - 16(1, 2)(1, 3)(4, 5) 

+ 10(1,2)(3,4)(5,6) . (3.30) 

Summing the three contributions we have 



~A 2 (1,2) = + 2(1,2) 3 -24(1,2) 2 (1,3)-8(1,2)(2,3)(3,1) 

+ 18(1, 2) 2 (3, 4) - 144(1, 2)(1, 3)(4, 5) + 96(1, 2) (2, 3)(3, 4) 

+ 60(1,2)(3,4)(5,6) + 24(1,2)(1,3)(1,4). (3.31) 



The left hand side of (|3~6| ) is: 



V[(l,2)] = (1,2)(1)(2)-4(1,2)(1)(3) + (1,2)(3)(4) + (1,2)(1) 2 -(1,2)(3) 2 . (3.32) 
Applying the 5 2 to the each term we get: 



^ 2 [(1,2)(1)(2)] = + (1,2)(2) 3 (1) + (1,2)(2) 2 (1) 2 -4(1,2)(1)(2) 2 (3) 

- (1,2)(1)(2)(3) 2 + 3(1,2)(1)(2)(3)(4) (3.33) 
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^ 2 [(1,2)(1)(3)] = + 3(1,2)(1) 2 (2)(3) + 2(1,2)(1)(2)(3) 2 -6(1,2)(1)(2)(3)(4) 

+ (1, 2)(1) 3 (3) + 2(1, 2)(1) 2 (3) 2 - 6(1, 2)(1) 2 (3)(4) 

- 9(1,2)(1)(2) 2 (3) + (1,2)(1)(3)(4)(5) (3.34) 

^ 2 [(1,2)(3)(4)] = + 2(1,2)(1) 2 (3)(4) + 2(1,2)(1)(2)(3)(4) + 8(1,2)(1)(3) 2 (4) 

- 16(1, 2)(1)(3)(4)(5) + 2(1, 2)(3) 3 (4) + 2(1, 2)(3) 2 (4) 2 

- 20(1, 2)(3) 2 (4)(5) + 20(1, 2) (3) (4) (5) (6) (3.35) 

^ 2 [(1,2)(1) 2 ] = + (1,2)(1) 2 (2) 2 + 2(1,2)(1)(2) 3 -4(1,2)(1)(2) 2 (3) + (1,2)(1) 4 
- 4(1, 2)(1) 3 (3) - 2(1, 2)(2) 2 (3) 2 + 6(1, 2)(1) 2 (3)(4) (3.36) 

^ 2 [(1,2)(3) 2 ] = + 2(1,2)(2) 2 (3) 2 + 2(1,2)(1)(2)(3) 2 + 4(1,2)(1)(3) 3 
- 12(1, 2)(1)(3) 2 (4) + (1, 2)(3) 4 - 6(1, 2)(3) 3 (4) 

3(1, 2)(3) 2 (4) 2 + 12(1, 2)(3) 2 (4)(5) . (3.37) 

We remind that at the fourth order the Wick rule produces 4 different topological 
graphs. Considering the invariance under permutation of the graphs labels we obtain: 

C[(l)(2)(3)(4)] = 3(1,2)(3,4), (3.38) 

C[(l) 2 (2)(3)] = 2(1,2)(2,3) + (1,2), (3.39) 

C[(l) 2 (2) 2 ] = 2(1,2) 2 + 1, (3.40) 

C[(l) 3 (2)] = 3(1,2). (3.41) 

We may now sum up the 5 contributions and apply the wick contraction C to the result. 
We obtain: 



-C[5 4 (l,2)] 



3 [ 2(1, 2) 3 - 24(1, 2) 2 (1, 3) - 8(1, 2)(2, 3)(3, 1) + 18(1, 2) 2 (3, 4) 
144(1, 2)(1, 3)(4, 5) + 96(1, 2)(2, 3)(3, 4) + 60(1, 2)(3, 4)(5, 6) 
+ 24(1,2)(1,3)(1,4)] 

= 3^[A 2 (1,2)], (3.42) 



which proves the theorem (|3.6j ) when the two operators in C5 4 and A 2 are applied 
to the graph (1,2). The extension to the general case is straightforward and already 
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included in the well definiteness of the operators S and A for general graphs as it was 
shown in section 6 of [|AC| . The (|3.6|) immediately imply the Q2.16j ). 
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